Driven Vortex States and Relaxation in Single Crystal Y Ba2Cu±0% 
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The vortex response to various ac-drives has been studied in high quality Y Ba2CuiOs, single 
crystals. Well within the vortex solid phase a re-entrant resistive state has been observed that is 
characterized by long relaxation times. Pulse current response and IV characteristics reveal that 
in this state the driven vortex system relaxes in a way to increase t he effective pinning force. A 
phenomenological model accounts for the observed features, in particular the clock-wise hysteresis 
of the IV-curves and the increase of the apparent critical current with applied current amplitude. 

PACS numbers: 74.60.Ge 



The behavior of the vortex system in type- II supercon- 
ductors in the presence of an applied drive current has 
been the subject of intense recent research. Of partic- 
ular importance for technical applications of supercon- 
ductors as well as for their fundamental understanding 
is the mechanism of depinning at a critical drive current 
which is believed to involve tearing of the vortex lattice 
and plastic flow 0,0 . Strong history effects and cur- 
rent induced transitions between high and weak pinning 
metastablc states f| have been reported. New relax- 
ation phenomena have recently been observed in NbSe2 
at sub-critical alternating pulse drives 0] a full theoreti- 
cal understanding of which has not yet emerged. 

Here we present magneto-transport measurements on 
single crystal Y BaiCu^O^ using various pulse and ac- 
drives. In the vortex solid state we observe a novel re- 
entrant, strongly non- linear response that is character- 
ized by long relaxation times. It is induced by ac-drives 
with amplitudes well below the dc-critical current. In 
contrast to the usually observed current-induced decrease 
in the pinning force, in this region of the phase diagram 
the driven vortex system relaxes in a way to increase 
the effective pinning force. A simple phenomenological 
model accounts for the principal observed features, in 
particular the clock-wise hysteresis of the IV-curves and 
the increase of the apparent critical current with applied 
current amplitude. We suggest that current-induced vor- 
tex entanglement is a possible mechanism underlying this 
behavior. 

The Y Ba2CuiOs, crystal used in this study has the 
dimensions a x b x c = 0.26 x 0.85 x 0.038 mm 3 and was 
grown in a self-flux method under high oxygen pressure, 
T c = 80 K (D . The occurrence of the various liquid and 
solid vortex phases was inferred from the temperature 
dependence of the resistivity for fields applied along the 
c-axis as shown in Fig. 1. Similar results have recently 
been reported by Qiu et al. Q. At a resistivity level of 
roughly 15% of the normal state resistivity at T c a sud- 



den drop is observed which smears out in high fields and 
that coincides with the onset of non-Ohmic transport as 
shown in the bottom inset of Fig. 1. This behavior is 
reminiscent to that found for un-twinned Y BaiCu-sOi-a, 
0, and in analogy we identify the onset of non-Ohmic 
behavior with the freezing/melting transition. The melt- 
ing line deduced from these data is well described by 
H m (T) — 32(1 — T/T C ) 1A and shown in comparison with 
the melting line of Y ' Ba-iCu-iOi _d || in the top inset 
of Fig. 1. Also shown is the scaling H m (T) / H cr with 
H cr = $o/(7s) 2 and 7 the anisotropy parameter and s 
the CwCVlayer spacing, which is expected for Josephson- 
coupled vortices fl . The melting lines of both materials 
can be approximately collapsed onto one curve. The val- 
ues of 7 = 7.5 and 12.5 have been used for Y BaiCu^O^-s, 
§ and YBa 2 Cu 4 8 @. 

Fig. 2 shows the field dependence of the resistivity 
at 70 K near the tail of the melting transition measured 
with a sinusoidal current with amplitude of 20 mA at 
various frequencies. The data are taken for increasing 
field after the sample has been zero-field cooled. Pro- 
nounced hysteretic effects between field cooled and zero- 
field cooled measurements will be described elsewhere 
|]lT|| . Well within the vortex solid state we observe a 
strongly frequency dependent re-entrant resistive state. 
The re-entrant behavior is reminiscent of the peak-effect 
phenomenon (T^l , and previous magnetic measurements 
|]lO| , [l3|Jl4[ at temperatures below 60 K have shown pro- 
nounced peak-effects superimposed on weak bulk pin- 
ning. However, the state observed here appears only for 
ac-drives: with dc-currents up to 60 mA (the highest cur- 
rent that can be passed without excessive contact heat- 
ing) no voltage response was detected. The frequency de- 
pendence of the resistivity at fixed fields is shown in the 
inset of Fig. 2. For fields in the re-entrant region a well 
defined cross-over from a strong increase roughly propor- 
tional to w 3 / 4 at frequencies below 500 Hz to a weaker 
w 1 / 4 or possibly logarithmic dependence at frequencies 
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above 500 Hz is observed. This defines a characteristic 
time scale of about 2 msec. In contrast, the melting tran- 
sition itself is essentially frequency independent in the 
frequency range studied here in agreement with results 
on NbSe 2 [E| and YBa 2 Cu 3 7 ^ d @. The long time 
constants in the re-entrant region indicate complex dy- 
namics possibly involving glassy relaxation and/or plastic 
deformations of the vortex system. 

The response to current pulses in the form of a sym- 
metric square wave is shown in Fig. 3. These data as 
well as those in Fig. 4 are obtained from an average over 
a large number of cycles using a digital storage oscillo- 
scope with 7 kHz bandwidth. They therefore represent 
the periodic response, the initial response on the first 
few cycles cannot be resolved. In the tail of the melting 
transition at 1.2 T the expected instantaneous voltage 
response in the form of a symmetric square wave is ob- 
served whereas in the re-entrant resistive region at 0.4 
T the voltage, after an instantaneous response on each 
current reversal, decays rapidly to essentially zero. The 
recurrence of the voltage response on every cycle shows 
that these data do not represent a current induced tran- 
sition from a weak pinning metastable state into a high 
pinning state as discussed in Ref. S. The data at 92 
Hz and 230 Hz indicate that the longer the time interval 
of constant current the further the voltage decays. The 
characteristic decay time of about 2 msec is consistent 
with the data in Fig. 2, and is independent of the ap- 
plied current amplitude ranging from 15 to 45 mA. This 
time dependence accounts for the strong frequency de- 
pendence shown in Fig. 2 and the absence of a response 
under dc-drives. Assuming that the data shown in Fig. 
3 represents the motion of the entire vortex system (this 
excludes channel flow) then an average travel distance of 
0.7 lira, can be estimated which corresponds to about 10 
vortex lattice constants. Very similar symmetric pulse 
responses have recently been reported for NbSe2 m- An 
oscillatory vortex response with periods much larger than 
those of the drive current has recently been observed in 
Y BaiCuzOi-d below the melting line |^7| . In contrast, 
the vortex response shown in Figs. 3 and 4 occurs on the 
time scale of the drive and represents therefore a different 
dynamic state. 

The IV-characteristics of Y Ba-JJu^O^ taken near the 
maximum of the re-entrant region with sinusoidal cur- 
rents are shown in the main panels of Figs 4a and 4b for 
70 K and 75 K. The highly non-linear voltage response 
(see right half of Fig. 4a) is characterized by a sizable 
hysteretic (that is, out-of-phase) component which has 
unusual clock-wise rotation. On the left half of Fig. 4a 
the in-phase component, V\ 1 and the out-of-phase com- 
ponent, V2, are shown separately. At current amplitudes 
around 15 mA the voltage response rises quickly. How- 
ever, with further increasing current amplitude the volt- 
age response does not follow the initial IV-characteristics 
but a characteristic that is progressively parallel shifted 



to higher current values: the apparent critical current in- 
creases with increasing amplitude of the applied current. 
At 70 K this unusual behavior persists to the highest 
attainable amplitudes. The progression of the in-phase 
component, Vi, of the IV-curves taken at 75 K is shown 
in Fig. 4b. The out-of-phase component at this temper- 
ature is strongly reduced as compared to the data in Fig. 
4a. At 75 K the usual IV-behavior is recovered for cur- 
rent amplitudes exceeding 40 mA. The apparent critical 
current as determined from the evolution of V\ at 70 K 
and 75 K and using a voltage criterion of 50 nV is shown 
in the lower inset of Fig. 4b. Except for current am- 
plitudes close to the threshold value the apparent critical 
currents for increasing and decreasing amplitude coincide 
and vary almost linearly with applied current amplitude. 

The response of the vortex system to an external elec- 
tromagnetic perturbation can be described on a macro- 
scopic level by the non-linear diffusion equation for the 
local current density |^8|: p,odi(r,t)/dt = \7 2 (p(j)j(r,t)). 
Numerical integration | 19| of this equation for current 
flow in a slap with a thickness much less than its width 
and non-linear resistivities of the form p(j) = po(j/j c ) n 
as well as p(j) = poexp(—j c /j) indicates that the results 
presented here can not be accounted for through current 
redistributions due to skin effect or inhomogenous criti- 
cal states. The data shown in Figs. 3 and 4 represent a 
novel dynamic state in which the vortex system adjusts 
itself under the influence of a driving force in a way to 
enhance the effective pinning force. We model this relax- 
ation behavior through a phenomenological extension of 
the hysteretic elastic creep theory ^0|. Within the elas- 
tic creep region the dominant current dependence in the 
IV characteristic is given by E = pj exp(—jth/\j\) where 
jth — Ucjc/ksT. Here j c is the critical current density, 
U c is the typical pinning energy, and p is the resistiv- 
ity scale. We adopt a value one for the glassy exponent 
and model the critical current relaxation by the simplest 
"equation of motion" for j th : dj t h/dt = -rjj th + g(\j\), 
where the constant 77 characterizes the rate of relaxation 
of jth to its equilibrium value and the function g de- 
scribes the effect of the applied current. In leading order 
we can expand g = go + S\j\+... For a sinusoidal varia- 
tion j — jo sin(wt), jth{t) contains aperiodic contribution 
and, in addition to gojr\ a constant off-set that is propor- 
tional to 5 jo. The parameters p,rj,go and S were deter- 
mined from fitting the data in Figs. 3 and 4 resulting in 77 
= 400 Hz, S = 8000 Hz, g Q = and p = 4.7 mfl cm. The 
fits are shown in the figures. The clock-wise hysteresis, 
the parallel shift of the IV-curves with current amplitude 
and the voltage decay for square-pulse currents, can be 
reproduced. In this model the clock-wise hysteresis arises 
because the periodic part of jth (t) reaches its maximum 
value with a delay with respect to the extrema of j(t). 
This delay is determined by the relaxation parameter 77 
and causes a smaller voltage on the decreasing branch 
of \j(t)\ as compared to the increasing branch. Similarly, 
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the shift of the IV-curves is caused by the constant off-set 
in jthit) that increases linearly with applied current am- 
plitude. This linear dependence is in fact observed in the 
apparent critical current (lower inset, Fig. 4b). We note 
that the apparent critical current extrapolates essentially 
to zero, consistent with go = 0. The value of rj = 400 Hz 
sets the time scale of the vortex response and is in good 
agreement with the frequency dependence shown in Fig. 
2 and the pulse-current response (see Fig. 3). However, 
the rapid voltage decay at short times (that is, high fre- 
quencies) in the pulse response as well as the threshold 
behaviors at the onset of the unusual IV-behavior and at 
the recovery of the usual behavior are not captured. This 
discrepancy may arise from the assumptions underlying 
this model, in particular, the linear expansion of <?(|j|), 
and a more complicated form of this function is required. 

We suggest that vortex entanglement |||] is a possible 
microscopic mechanism underlying the observed behav- 
ior. An applied drive current can cause the forced entan- 
glement of vortex lines which will lead to a substantial 
suppression of vortex creep in much the same way as dur- 
ing work hardening the increasing number of dislocations 
effectively leads to a blockage of their motion. This argu- 
ment is in qualitative agreement with the simulations of 
vortex motion in three-dimensional layered superconduc- 
tors modeled by the time-dependent Ginzburg-Landau 
equation. As a measure of entanglement we calculated 
the diffusion of vortex lines under the influence of a drive 
current, (\r(z) — r(z')\ 2 ) = D\z — z'\. Here, r(z) is the 
position of the vortex line, 2 the coordinate along the field 
direction and D the "diffusion" coefficient. Our simula- 
tions indicate a sharp increase of the entanglement near 
the onset of vortex motion, see top inset of Fig. 4b. 

In conclusion, the states of vortex dynamics in 
Y BaiCu^O^ have been investigated using magneto- 
transport measurements. The melting line H m (T) sepa- 
rating the vortex liquid and vortex solid states has been 
determined from the onset of non-Ohmic behavior. H m 
shows good scaling according to H m ~ 7~ 2 with the 
results of YBa2CusOj-d- Well within the solid state 
a re-entrant strongly current and frequency dependent 
resistive state is observed. Pulse-current response and 
IV-characteristics reveal that in this region of the phase 
diagram the driven vortex system relaxes in a way to in- 
crease the effective pinning force. A simple phenomeno- 
logical model accounts for the observed clock-wise hys- 
teresis of the IV-curves and the increase of the appar- 
ent critical current with applied current amplitude. We 
argue that forced vortex entanglement due to the drive 
current is a possible microscopic mechanism underlying 
this dynamic vortex state. 
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T(K) 

FIG. 1. Temperature dependence of the resistivity of 
Y Ba2CuiOg, in various fields applied parallel to the c-axis. 
Bottom inset: Onset of non-Ohmic behavior at the melting 
transition. Top inset: Comparison of the melting lines H m (T) 
(solid symbols) and scaled melting lines H m /H cr (open sym- 
bols) of Y Ba2Cu^.O-7-d and Y Ba-zCu^O%. 
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FIG. 2. Field dependence of the in-phase resistivity at 
70 K measured with ac-currents of various frequencies. Inset: 
Frequency dependence of the resistivity at two fields in the 
re-entrant region and in the tail of the melting transition. 
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FIG. 3. Voltage response to currents in the form of 
symmetric square-pulses at 70 K for magnetic fields in the 
re-entrant region and in the tail of the melting transition. 
The heavy broken line is a fit according to the phenomeno- 
logical model. 
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FIG. 4. (a) Current-voltage characteristics at 70 K and 
0.4 T for various amplitudes of the applied sinusoidal current. 
The right half of the figure shows the total voltage response. 
The clock-wise sense of rotation is indicated by the arrows. 
In the left of the figure the in-phase and out-of-phase compo- 
nent are shown separately. The lines are fits according to the 
phenomenological model; (b) Current-voltage characteristics 
for the in-phase component Vi at 75 K and 0.25 T. Bottom 
inset: Amplitude dependence of the apparent critical current 
at 70 K and 75 K as determined from a 50 nV criterion. Top 
inset: Current dependence of the vortex diffusion coefficient. 



